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Classical variational principles can be deduced from quantum variational principles via formal 
reparameterization of the latter. It is shown that such reparameterization is possible without in¬ 
voking any assumptions other than classicality and without appealing to dynamical equations. As 
examples, first principle variational formulations of classical point-particle and cold-fluid motion are 
derived from their quantum counterparts for Schrodinger, Pauli, and Klein-Gordon particles. 


I. INTRODUCTION 

It is commonly known that variational methods are a 
powerful tool for studying the dynamics of various phys¬ 
ical systems, e.g^ quantum molecular dynamics [ll, 
fluid mechanics [3, jj, classical motion of single parti¬ 
cles and collective processes in plasmas ii , and wave 
propagation in both linear and nonlinear media [M3|. 
Variational formulations are advantageous as they lead 
to dynamical equations in a manifestly conservative form 
derived from a single scalar function, a Lagrangian (or 
Lagrangian density). The fact that these equations can 
be approximated robustly and self-consistently by ap¬ 
proximating just one function makes the method particu- 
larly attractive for reduced calculations [1, [Ill, [T^ • How¬ 
ever, so far, exact Lagrangians have been obtained largely 
heuristically or ad hoc (l3l| . Such approaches tend to ob¬ 
scure the physical meaning of the results and limit their 
applicability, while regular ways to deduce Lagrangians 
rigorously from first principles are yet to be found. 

Here we show that classical variational principles 
(VPs) can be deduced from quantum VPs, which are 
well known, via formal reparameterization of the lat¬ 
ter. Such reparameterization is possible without appeal¬ 
ing to dynamical equations and without invoking any 
assumptions other than classicality. This distinguishes 
our theory from the existing variational formulations of 
the quantum-classical correspondence, which are more 
restrictive H- Also as a complement to those formula¬ 
tions, we consider both single-particle and fluid VPs and 
rigorously explain how they are connected through quan¬ 
tum VPs. Classical-fluid Lagrangians flow as the semi- 
classical limit of the fundamental quantum Lagrangian 
(FQL), Eq. ([3]), and the point-particle Lagrangians are 
then yielded as corollaries for narrow but otherwise ar¬ 
bitrary wave packets. This approach enables the first 
principle classical Lagrangian description of most gen¬ 
eral, relativistic vector particles (e.g., a Dirac electron) 
and, similarly, the geometrical optics (GO) description 
of any vector waves, as we discuss in a companion paper 
[T^. The present paper is intended as an introduction 
to such calculations. Thus, below, we primarily focus on 
systematization of VPs for commonly known systems, in¬ 
cluding Schrodinger, Pauli, and Klein-Gordon particles. 

We explicitly show how to deduce classical Lagrangians 
for these particles and the corresponding fluids from 
the FQL. In particular, we show that the expression 


for the Bohm quanturn potential dill , the so-called 
Weizsacker correction [l8|, the Madelung equations [l^ . 
and the Ghen-Sudan Lagrangian density [20j all emerge 
naturally within our unifying theory as special cases. We 
also obtain an alternative, manifestly Lagrangian repre¬ 
sentation of Takabayasi equations for a Pauli particle. 
In contrast to the original theory (2lj |. our model yields 
the dynamics of the spin vector S directly from the VP 
and employs two, rather than three, equations; hence, 
S' = 1/2 is ensured irrespective of initial conditions. 

The paper is organized as follows. In Sec. m we define 
the basic notation. In Sec. nil we introduce the general 
formulas. In Sec. IlYl we discuss a Schrodinger parti¬ 
cle and the Madelung equations. In Sec. El we discuss 
a Pauli particle and the variational formulation of Tak¬ 
abayasi equations. In Sec. ED we discuss a Klein-Gordon 
particle and rederive the Ghen-Sudan fluid Lagrangian 
density. In Sec. IVHl we summarize our main results and 
outline future applications of the proposed theory. 


II. NOTATION 


The following notation is used throughout the paper. 
The symbol “=” denotes definitions. We use natural 
units, so the speed of light equals one (c = 1), and so 
is the Planck constant {h = 1). The Minkowski met¬ 
ric is adopted with signature (-,-!-,+,+), so, in partic¬ 
ular, d'^x = dtd^x. Generalizations to curved metrics 
are straightforward to apply [^ . Greek indexes span 
from 0 to 3 and refer to spacetime coordinates, x^, with 
corresponding to the time variable, t\ in particular, 
dfi = djdx^. Latin indexes span from 1 to 3 and denote 
the spatial variables, x® (except where specified other¬ 
wise); in particular, di = djdx^. Summation over re¬ 
peated indexes is assumed. Also, in the Euler-Lagrange 
equations (ELEs), the notation “da denotes, as usual, 
that the corresponding equation was obtained by extrem- 
izing the action integral with respect to a. 
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III. BASIC EQUATIONS 
A. Fundamental quantum Lagrangian 


The dynamics of a quantum particle is governed by the 
least action principle, 


Oq 

> 

II 

0 

(1) 

where A is the action integral. 


II 

<1 

(2) 

and £ is the Lagrangian density. 

Like for any other 


nondissipative linear wave (2^ . the Lagrangian density 
of a quantum particle can be expressed as [23j, 

£= (3) 

which is what we term (the density of) the FQL. Here 
H is some Hermitian operator called Hamiltonian, ?/; is 
a complex vector field (“state function”), and ijj^ is its 
adjoint. Treating ip and ip^ as independent functions [^ . 
one obtains from Eq. © two mutually adjoint ELEs, 


which one may recognize as Hayes’s representation of the 
Lagrangian density of a GO wave [2^. Treating 6 and 
I as independent functions, one hence arrives at the fol¬ 
lowing ELEs: 

S9: dtl + V ■ (Iv) = 0, (9) 

SI: dt0 + H{t,x,V0) =0, (10) 

where v = dpH{t, x, p) is the group velocity evaluated at 
the wave vector p(t, x) = V0. Equation ([9]) is a continu¬ 
ity equation known as the action conservation theorem 
(ACT). Equation (flUl) is the Hamilton-Jacobi represen¬ 
tation of the local dispersion relation, 

f = iJ(t,x,p), (11) 

where £ = —dt0 is the local frequency. 

To the extent that a wave packet is well localized such 
that it is meaningful to describe its dynamics as the dy¬ 
namics of the packet’s geometrical center, the continuous- 
wave description can be replaced with a simpler, point- 
particle model. In this case, one can approximate the 
action density with a delta function, 

I(t,x)=<5(x-X(<)), (12) 


Sip^ : idtip = Hip, (4) 

Sip : —idtip^ = {Hip)"^. (5) 

It is also possible to derive other, different equations 
from Eqs. © if £ is reparameterized, i.e., if other inde¬ 
pendent functions are chosen instead of ip and ip"'^. This 
can be done in general, using the fact that the Hamil¬ 
tonian always can be expressed as H = H(t, x, p) (2^ 
(here H is some matrix function, x = x is the position 
operator, p = —fV is the momentum operator, and the 
standard coordinate representation is assumed), so 


so the LD can be replaced with just a point-particle La¬ 
grangian, L = J £d^a:. This is done as follows. Let us 
express the complete action integral, A, in the form 


A = — JX(t, x) dt0{t, x) 


— y I(t, x) iJ(t, X, V0(t, x)) d^x . (13) 


A2 


-2 = ^ [p^\dtip) - {dtip'^)ip] - ip^H{t,yL,-iV)ip. (6) 

Specifically, the classical reparameterization of this La¬ 
grangian density is performed as follows. 


Then, substituting Eq. m leads to 
Ai = — y S{x — X(t)) dt0{t, x) dP^x 

= J [clti5(x - X(t))] 0{t, x) d'^x 


B. Single-particle Lagrangian 

First, let us consider the simplest classical reparam¬ 
eterization of Eq. For a scalar particle, the state 

function can be expressed as 


Ip = ae 


i0 


( 7 ) 


where the amplitude a = 'Jl and the phase 0 are real 
scalar functions of {t, x). Let us assume the semiclassical 
limit, i.e., that X(t,x) is slow compared to 0(t,x). Then, 
Eq. m gives m 

( 8 ) 


= J [-V5(x-X(t))] •X(t)6»(t,x)d^x 
= / - J 0{t,x)VS{x-X{t))d^x ■X{t)dt 
'' S{K-X{t))V0{t,x.)d^x ■X{t)dt 
= JP{t)-X{t)dt, (14) 

where P(t) = V6>(t,X(t)). Similarly, 

A. = -/ 


-I[dt0 + H{t,x,V0)], 


H{t,X{t),P{t))dt. 


(15) 
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Thus, the total action integral is expressed as A = J Ldt. 
The function L serves as a point-particle Lagrangian and 
is given by 

L = P-X-i/(t,X,P). (16) 

One may recognize this L as the standard classical La¬ 
grangian, where X and P serve as a canonical coordinate 
and a canonical momentum, respectively. Treating them 
as independent variables then leads [1^ Sec. 43] to ELEs 
matching Hamilton’s equations in the canonical form, 

5P: ±^dpH{t,X,P), (17) 

5X: P =-dxH{t,X,P). (18) 

(Hence, £ and p can be understood as local canonical 
energy and momentum.) Equations and m also 
serve as ray equations for Eqs. (ED and m M- 

It is to be noted that our derivation does not require 
interpreting ip as a probability and holds for any H. It 
is also to be noted that our point-particle equations are 
derived without assuming a specific profile of I{t, x). All 
we require instead is that this function be narrow enough 
(yet wide enough compared to the particle wavelength; 
otherwise, the GO approximation does not apply). 


£ given by Eq. EID can be viewed as the Lagrangian 
density of a cold classical fluid. (Specific examples will 
be discussed in the next sections.) It is to be noted, like 
we already did before, that interpreting ip as a probability 
amplitude is not needed within this approach. 

Thermal corrections for isentropic fluid also can be in¬ 
troduced readily, at least ad hoc. This is done by adding 
to £ a term —nU{t,x,n), which leads to 

£ = —n[dtd + + U(t,x,n)]. (24) 

Erom comparing Eqs. EID and (l24l) . it is seen that U 
has the meaning of the per-particle internal energy; then 
n^dnU serves as the pressure 

Outside the ad hoc approach, there exists other vari¬ 
ational methods to obtain corrections to the cold-fluid 
variational principle that we derived. In Ref. [2^, the 
Thomas-Eermi energy functional is minimized in order to 
obtain the quantum hydrodynamic equations with ther¬ 
mal gradient corrections included. Also, the L^rangian 
density for more general, non-isentropic fluids [3| can be 
deduced from the FQL. Details will be reported in a sep¬ 
arate paper. Below, we will completely ignore thermal 
effects for clarity. 


C. Fluid Lagrangian density 


IV. SCHRODINGER PARTICLE 


Now consider an ensemble of particles described by 
{9q,Iq), q = Assuming the energy of parti¬ 

cle interactions (except through average fields included 
m H) is negligible, the ensemble Lagrangian density £ 
equals the sum of single-particle Lagrangians 2q given 
by Eq. (ED- Let us also assume that particles are cold. 
This means that all are about the same, and thus so 
are all £q\ i.e.. 


V9q ^ V9, dt9q SS dt9, (19) 

where 9 is some average phase. Then, the ensemble is 
described by the same £ as the one given by Eq. ®, 
except I must be replaced with the particle density. 


Let us specifically discuss a Schrodinger particle, i.e., 
a nonrelativistic charged particle governed by 

^ (-iV - qAf + q(p. (25) 

Here m and q are the particle mass and charge, A = 
A(t,x) is the vector potential, and (p = is the 

scalar potential. Substitution of Eq. 0 leads to the fol¬ 
lowing exact Lagrangian density. 


£ = -I 


dt9 +—{V9 - qAf + qp, 
Zm 


lEs = 


(Va)2 (VI)2 


2m 


8ml 


(26) 

(27) 


N 

n = Y^Xq. ( 20 ) 

In other words, the ensemble’s Lagrangian density is 
given by 


£ = —n[9t0-I-iL(<, X, V0)]. (21) 

The corresponding ELEs are 

59 : dtn + V • (nv) = 0, (22) 

5n-. 9(0-I-iL(t, X, V0) = 0, (23) 

and V serves as the flow velocity. These equations provide 
a complete description of the ensemble as a cold fluid, so 


where Wg can be recognized as the Weizsacker correction 
M- The corresponding ELEs are 

69: dtl+V ■ (Iv) = 0, (28) 

5X : dt9 t;— (V 9 — qA)^ q^p Q = 0, (29) 

Zm 

where v = (V0 — qA)/m, and 


Q = 


Zmy/X 


(30) 


is recognized as the Bohm quantum potential EIH, 
which is obtained when varying the J kEs d^a: term of the 
action. Equation (1^51) represents the ACT, and Eq. ESD 
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represents a quantum Hamilton-Jacobi equation. Taking 
the gradient of the latter readily yields 

m{dt + V • V)v = ( 7 (E -I- V X B) — VQ- (31) 

Equations (1281) and dSID combined together coincide 
with Madelung equations [l^. Hence, £ given by 
Eq. (|26]) can be identified as the Madelung Lagrangian 
density. The equation for cold classical fluid are obtained 
by replacing I with n, as in Sec. IHI Cl and also by ne¬ 
glecting Q in Eq. (1^ or, equivalently, ITs in Eq. (l26ll 
(which constitutes the semiclassical, or GO approxima¬ 
tion). Hence, 


£ = -r 


dtO + - qAf 

2m 


q(p 


(32) 


can be identified as a Lagrangian density of a cold clas¬ 
sical fluid, in agreement with Ref. Q. According to 
Sec. IHI B1 the classical point-particle motion is then gov¬ 
erned by Eqs. (13 and (flSl) with the Hamiltonian 

H{t, X, P) = ^ [P - qA{t, X)]2 + q^{t, X). (33) 

Zm 


V. PAULI PARTICLE 

Now let us consider a Pauli particle^ _ i.e., a non- 

relativistic spin-1/2 particle governed by 


H = ^(-*V - qAf + qq,-^cT-B, 

where cr = (crx,cry,az) are the Pauli matrices, 
'0 1\ /O -i\ fl 0 


1 0 


i 0 


0 -1 


(34) 


, (35) 


and B = V X A is magnetic field. [Replacing the Bohr 
magneton qjiZm) in Eq. (IMl) with a different constant 
would not significantly affect the below calculation.] The 
state function ip has two complex components; i.e., 

where ai _2 and 9i 2 are real. Thus, it can be cast in the 
following general form, p; = QVx, with 

/ n cos(C/2)\ 

Here 9 ^ {9i+ 02 ) 72 , d = 92 —9i, I = V'V, C is real and 
defined via 


oi = acos((//2), 02 = asin(C/2), 


(36) 


and a = \/T. Hence, the Lagrangian density can be 
expressed as 


£ = y Widt-q) - {dtq^)q\ - I{dt9) - Wp 


-I 


2m 


(V0 - qA)"^ + q‘p + UsG 


(37) 


where d* = -I- v • V is the convective derivative associ¬ 

ated with velocity field v = (V0 — qA)/m, and 


UsG = B)77, 


(38) 


where we made use of the fact that rj^rj = 1. The term 
UsG is the Stern-Gerlach energy. It can be expressed 
conveniently as Usg = —{<11^^) S • B, where 


^ ^ /sin (/ cos d' 

S = - rj^crr] = — sinC sind 
2 


(40) 


is the spin vector, S = 1/2. Expressing also the other 
terms through the four independent variables, {9,1, d, /), 
one gets 


£= -I 


dtO + 7^ (V0 - qAf 

2m 


qif 


- ^{dtd) cosC- — S{d,() ■ B 

2 m 


-VPp, 


Wp = 


(V^)^ 

8toI 


— [(vc)^-f(Vd)^]. 


(41) 


Hence, the following four ELEs are yielded. The first 
one is the ACT, 

69: dtl+V- (IV) = 0, (42) 

where the flow velocity is given by V = v -|- u, and 

u = - (Vryi)??] = (Vd)cosC. 

The second ELE is a Hamilton-Jacobi equation, 

SI : dt9+^ (V0 - qAf +qip + Q-^S-B 

2m m 

= i(dtd)cosC-^[(VC)2 + (Vd)2]. (43) 

Finally, the remaining two ELEs are 

Sd : 9t(Icos/) + V • (Ivcos/) 

= 7^ V • (IVd) + ^ (a,,S) • B, (44) 
2m m 

and 

JC : I {dtd) sin C = 7 ^ V • (IV/) + — (9^8) • B. 

2m m 

(45) 

As shown in appendix, Eqs. (H51) - (H5l) rep resent an al¬ 
ternative form of Takabayasi equations [2l| , 

m{dt + V • V)V = q(E -P V X B) - VQ 

+ ^ ^ dj [I(V5,)(9^^*)] , (46) 
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^tS + CV■V)S = —SxB + ^Sx^^(I^^S) (47) 
m mX 

(which were also studied recently in Refs, How¬ 

ever, our equations have a number of advantages: (i) they 
have a manifestly Lagrangian form; (ii) the spin dynam¬ 
ics is described by two rather than three equations; (iii) 
S' = 1/2 is ensured irrespective of initial conditions. 

The Lagrangian density for cold fluid is obtained by 
replacing X with n, as in Sec. Imcl and by assuming 
that the real angles, which describe the particle spin, are 
about the same; i.e., 'dq ~ d and Cg ~ Cj where "d and C 
are the corresponding average angles. Hence, 


representation and also a detailed comparison with re¬ 
lated point-particle models of spin electrons is discussed 
in a companion paper , where our results are extended 
to relativistic particles. 

VI. KLEIN-GORDON PARTICLE 

It is instructive to discuss also Lagrangian densities of 
other types. For example, let us consider 

= (56) 


£ = -r 


dte 


-(Vd-,A)^. 


qip 


- (dtd) cos C - — S(d, C) • B 

z m 


- (48) 


which corresponds to a Klein-Gordon particle [i^ 
i.e., a relativistic spinless particle governed by 

+ m'^)^ = 0. (57) 


can be identified as a Lagrangian density of a cold fluid 
comprised of Pauli particles. The point-particle La¬ 
grangian is derived much like in Sec. IHIBI We neglect 
the Bohm quantum potential and also assume that "d and 
C are approximately homogeneous in the regions in which 
I(t,x) is nonzero. Specifically, one obtains 


Here — qA^^, = A^{x'') is an electro¬ 

magnetic four-vector potential (so (p = —Aq), and ip is 
complex scalar function. [The factor (2 to)“^ in Eq. ((56ll 
is added to ensure that ip is normalized such that ip*ip 
relates to the action density most transparently; see be¬ 
low.] Let us represent the state function as 


L = P-X-bidcosC-Z7(t,X,P,d,C), (49) 

where the Hamiltonian H is given by 


= X (P - qA{t,X)f 
2m 

+ 5 <^(i,X)-^S(d,C)-B(f,X). (50) 

m 

The corresponding ELEs are 

dP: X = apiZ(t,X,P,dl,C), (51) 

dX: P = -5xZ7(LX,P,d,C), (52) 

dd : sinC = —{2q/m) (9^S) • B, (53) 

d(/: d sin(/=-|-(2g/m) (9(;S) • B. (54) 

These equations also yield the following equation for S, 
S =-^S X B(t,X), (55) 


iP = e^^s/To (58) 
and adopt the semiclassical limit, like in Sec. mu Then, 

+ -{Ve-qKf-m% (59) 

This has the form of Whitham’s Lagrangian density m, 

£ = D(-5td, Vd;t,x)Io, (60) 

e p 

in which 

X = Xode1){£,p]t,x) =jXo (61) 

serves as the action density, and 7 = (£ — q(p)/m. This 
yields the ACT and a Hamilton-Jacobi equation, 

dd : dtX+V- (Iv) = 0, (62) 

5X : {dtO + qp)'^ = (V9 — qA)"^ + m'^, (63) 


as can be checked by direct substitution. Also, even with¬ 
out such calculation, the fact that S is perpendicular to 
S is anticipated from the definition of the spin vector, 
Eq. (1401) . according to which S must remain constant. 

Notice that, within our approach, the spin preces¬ 
sion equation flows directly from the VP. This distin¬ 
guishes our approach from other calculations (e.g., in 
Refs, [131), where the equation for S is derived from 
a postulated Poisson bracket for the classical spin vari¬ 
ables. We also notice that our Lagrangian can be repre¬ 
sented in terms of complex variables, in a form akin to 
but different from those in Refs. [38l - l44| . This complex 


where the flow velocity v is given by 

V0-qA 

V = -. 

ni'y 


(64) 


Hence one may recognize 7 as the Lorentz factor [be¬ 
cause £ serves as the canonical energy (Sec. IHI|) and, 
from Eq. (IS51) . it is seen that £ — qp is the kinetic en¬ 
ergy] , so Xq is recognized as the proper action density. 

Erom Eq. (1551) . one gets 


-ftd = ±V(Vd-(?A )2 + w? + qip. (65) 
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To ensure that the previously adopted semiclassical ap¬ 
proximation is satisfied, one of the two solutions must 
be excluded. We will consider positive-energy particles, 
i.e., assume the plus sign. Then the wave is of a scalar 
type, and, as it is well known [2^, Whitham’s Lagrangian 
density can be rewritten in the Hayes’s form ([8]); namely, 

£ Ri -X{dt9 + a/[p - gA(t,x )]2 + 711 ^ + g(^(t,x)}. 

( 66 ) 

As in the previous sections, the Lagrangian density of a 
cold (yet now relativistic) classical fluid is obtained by re¬ 
placing X with n. The corresponding ELEs coincide with 
Eqs. (1^ and (1^ . and taking the gradient of the latter 
leads to the well known relativistic momentum equation, 

[dt -I- (v • V)](to 7 v) = g(E -|- v x B). (67) 

Also, the classical point-particle motion is automatically 
seen to be governed by Eqs. (HZD cinci (jlSp witli 

H{t, X, P) = V[P - 9A(t, X)]2 + q^{t, X). (68) 


Finally, let us discuss the interactions of particles with 
self-consistent fields. The Lagrangian density of such 
interactions is obtained, as usual, by adding to £ the 
Lagrangian density of the vacuum field [i^ Sec. 11.5], 


1 . VT u *n2 |Vx A|2 _ 

£EM = ^(-V7^-a*A) = 

where E and B are electric and magnetic fields. This 
gives 

£e = £em - '^ns[dt9s + L7s(i,x, Vds)], (69) 


where the sum is taken over all species. In particu¬ 
lar, substituting here the Hamiltonian (1681) automatically 
yields the well-known Chen-Sudan Lagrangian density of 
cold relativistic plasma [l^. It is to be noted that, in 
the form ((6^ . the obtained £ 1 ; can be treated as the La¬ 
grangian density of warm plasma too, if particles with 
different energies are treated as different species; for ex¬ 
ample, see Ref. [HI- 


We showed that, for instance, the expression for the 
Bohm quantum potential, the so-called Weizsacker cor¬ 
rection, the Madelung equations, and the Chen-Sudan 
Lagrangian density all emerge naturally within our uni¬ 
fying theory as special cases. We also obtained an al¬ 
ternative, manifestly Lagrangian representation of Tak- 
abayasi equations for a Pauli particle. Our model yields 
the dynamics of the spin vector S directly from the VP 
and employs two, rather than three, equations; hence, 
S' = 1/2 is ensured irrespective of initial conditions. 

Our results can be viewed as a generalization of and 
a complement to reparameterizations of the FQL that 
were reported in literature earlier [T3 |. In addition to 
this, the new approach has several important applica¬ 
tions. First of all, it enables a first principle classical 
Lagrangian description of most general, relativistic vec¬ 
tor particles such as a Dirac electron, which is a long¬ 
standing problem (see, e.g., Ref. 0). Second, the same 
theory is applicable as is to a GO description of clas¬ 
sical vector waves, since the fundamental Lagrangian of 
classical waves is identical to the FQL This allows 
accounting for polarization effects and mode conversion 
directly in ray equations. This also allows to extend the 
recent studies of ponderomotive forces on scalar waves 
and particles [l^ to more general, vector waves and par¬ 
ticles. Third, calculations of even purely classical dynam¬ 
ics can be simplified through the application of quantum 
Lagrangians. This is because, in contrast to classical 
equations for point particles, semiclassical equations are 
linear and thus sometimes are easier to work with. These 
and other applications of the theory presented here will 
be discussed in follow-up papers, including Ref. fisj ]. 

The authors thank J. W. Burby and N. J. Fisch for 
valuable discussions. The work was supported by the 
NNSA SSAA Program through DOE Research Grant 
No. DE274-FG52-08NA28553, by the U.S. DOE through 
Gontract No. DE-AG02-09GH11466, and by the U.S. 
DOD NDSEG Fellowship through Contract No. FA9550- 
ll-C-0028. 


Appendix A: Takabayasi equations 


VII. DISCUSSION 

We showed that classical VPs can be deduced from 
quantum VPs, which are well known, via formal repa¬ 
rameterization of the latter. Such reparameterization is 
possible without appealing to dynamical equations and 
without invoking any assumptions other than classical- 
ity. Classical-fluid Lagrangians flow as the semiclassi¬ 
cal limit of the FQL, and the point-particle Lagrangians 
are then yielded as corollaries for narrow but other¬ 
wise arbitrary wave packets. We explicitly performed 
these calculations for several commonly known systems, 
namely, Schrddinger, Pauli, and Klein-Gordon particles. 


In this appendix, we show that the ELEs derived in 
Sec. El for a Pauli particle agree with Takabayasi equa¬ 
tions (H51) and (HTl) . 


1. Momentum equation 


To derive the momentum equation, we begin byrewrit- 
ing the “spin stress” introduced by Takabayasi |2l| . 


n = -^d,[x{vs.){d^s^)], (Ai) 
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in Eq. (gSI). Specifically, by substituting Eq. go]) for S 
in Eq. dUD , one can write 

n = [I(VC)(5^C) +Isin2 C (Vi9)(a^-^)]. 

(A2) 


Also, let us use Eq. (l45l) to express (it'd and substitute 
the result in the right hand side of Eq. (H51) . Taking the 
gradient of the resulting equation yields 


J 


mdtv = 9(E + V X B) - VQ + ^ V [S • B + cot C (9^8) • B] 


cot C 
Ami 


V • (IVC) 


- ^ V [(VCf + (Vd)2] . (A3) 


Proceeding with the calculation of the conservation equation for tou, one obtains 


mdt\i = - ^9t[(Vd) cosC] 

= - i (Vi?) (cos C) - ^ cos C V {dt'd) 


= - (Vi?) 

— cos C V 


-i (v V) cos C + V • (lu) 


V • (IVi?) + ^ (d^S) ■ B 
4mx m 


_i (v V)i? +- - - 

2 ' Ami sin ( 


V • (ivc) + 


<7 


TO sin C 


(dcS )■B 


771U fVl?) 

= - to(v V)u - m(Vviy + — V • (lu) - V • (IVi?) 

I Ami 


--^[(9^S)-B] Vi?-cosCV 




TO sin C 


(d^S) ■ B 


— cos C V 


1 


4toI sin ( 


V • (ivc) 


COS C (^'i9) 

= - to(v V)u - to(u • V)v + gu X B + —X (Vi?)V • [(Vi?)IcosC] - V • (IVi?) 


- — [(9^,8) •B](Vi?) -cosCV 


TO sin ^ 


Ami 
(c?c8) • B 


4toI 


— cos C V 


4toI sin C 


V • (ivc) 


(A4) 


where we used Eqs. (HI) and (1151) . Let us add Eqs. (IA3I) and (IA4I) and also add (u • V)u to both sides in order to 
complete the convective derivative V = dt+'V • V. One then gets 


toPV =g(E + V X B) - VQ + — S\VB^) + to(u • V)u - V • [IVQ - ^ V[(VC)^] 

m AmX 8m 

^ (Vi?) V • (IVi?) - (Vi? • V)(Vi?) - (Vi?) [(Vi?) • (VC)] 


Ami 


Am 


2m 


= g(E + V X B) - VQ + ^ S\VB,) - d, [l[VC)[d^C) +Isin^ C (Vi?)(a^i?)] 

m 4:mZ 

= g(E + V X B) - VQ + ^ S\VB,) - -X [l(V5,)(9^.8*)] , 

m mX 


(A5) 


where we substituted Eq. (IA2I) at the end. Hence, we see that the ELEs derived in Sec. El lead to Eq. (1151) . 


2. Spin equation 

To derive the spin equation, let us start by re-expressing the “spin torque” introduced by Takabayasi 


M = ^[sxdj[id^s)], 


(A6) 



































in terms of 1 !) and C- Specifically, we get 


M, = ^ SyV ■ {TVS,) - ^ 5 ,v • {IVSy) 

sin^sin-d V • [2iV(cosC)] — ^ cos^ V • [IV(sin^sini?)] 


Ami 


Ami 


Ami 
sini? 


sin ^ sin-d V • [I(VC) sin^] — ^ cos^ V • [I(VC) cos C sin-d] — ^ cos^ V • [I(Vi?) sin ^ cos-d] 


Ami 


Am 


{V^Y sin C cos C - ^ V • (IVC) 


1 


— -— (VC)^ sin ^ cos C sin-d — -— cosC(VC) • V(cosCsinz?) — 


4m 
sini? 
4m 


4m 

{VdY sin C cos C - V • {IVC) 


Ami 


Ami 


cos C cos dV ■ [I( Vz?) sin 


Ami 


cot C cos-d V • [I(Vz?) sin^ C]) 


(A7) 


My= — S,V ■ {IV SY -—^S,V- {IVSY 
ml ml 

1 


1 

ml 
1 

4^1 ■ L- • V .j 4 ^ 2 : 

cosC V • [I(VC) cos C cos-d] — ^ cos^ V • [I(Vi?) sin ^ sin z?] + ^ sin ^ cos r? V • [I(VC) cos^] 


cosCV • [IV(sin(^cosi?)] — -—— sinf^cosz? V • [IV(cosC)] 


4ml 
cos I? 


4ml 


4m 


1 


{V-dY sin C cos C - :y V • {IVQ 


1 


Ami 


Ami 

cotCsinz?V • [I(V-d) sin^ (^], 


(AS) 


M, = ^ S^V ■ {IV Sy) -^SyV- {IV SY 
sin () cos "0 V ■ [IV(sin ^ sin i?)] — 


Ami 
sin ^ cos C 


4ml 


sin ^ sin z? V • [IV (cos C cos -d)] 


2m 


(VC) • (v^) 


sin C , sin C 


Ami 


(VI) • {VY + 


4m 


{V-d) 


1 


Ami 


V • [I(Vz?) sin^ C]- 


(A9) 


Hence, the conservation equations for the spin vector S can be expressed as follows: 


lASx = - i9t (sin C cos I?) + - (V • V)(sinCcosr?) 

= - cos C cos -d PC ~ 2 ^ ^ 

= — cot C cos d VSz — - sin -d sin C P'd 

= — cote cos d I -—— V • [I(Vd) sin^ C] + “ (5^S) • B 

4ml m 


— sin-d 


i sinC (u • V)t9 + ^ V • (IVC) + ^ (acS) • B 


= — cot C cos d 


V ■ [I(Vt9) Sin2 C] + ^ (a^S) • B 
4ml m 


— sind 
1 


(Vd)2 sin C cos C + ^ V • (IVC) + - (acS) • B 
4m 4ml m 


Ami 


cot Ccosd V • [I(V-d) sin^ C] + sind 


^(Vd)^sinCcosC-^V.(IVC) 


(S X B) 


= ^ [s X a,(ia^s)],, + ^ (s X B)„ 

mX m 


(AlO) 
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where we substituted Eq. (IX7l) : also, 


VSy = 7 : dt (sin C, sin 1 ?) + - (V • V) (sin C, sin -d) 


= - cos C sin d + 2 C cos d V-d 
= — cot C sin i? VSz + - cos d sin Q Vd 


= — cot C sin d 

'1 


V • [I{V^) sin^ C] + ^ (d^S) ■ B 
4rnX m 


COST? 

1 

AmX 


2 sinC (u • V)d + ^ V • (XVC) + ^ (9cS) • B 


cot (^sind V • [X(V'd) sin^ — cos'd 


^(V^fsinCcosC-^V.(XVC) 


= ;^[Sxa,(X9^S)], + ^(SxB)„ 


where we substituted Eq. (IA8E Einally, the equation for Sz is obtained directly from Eq. (l44l) . 


m 


(S X B) 


(All) 


T[dtSz + (v • V)5,] + Sz[dtT + V • (Xv)] = ^ V • (XVz?) + ^ (9^8) • B. 

4m m 


(A12) 


Using the ACT, one then gets 

VSz [T{n ■ V)Sz + • (Xu)] + -^ V • (XVd) + ^ (9^8) • B 

X 4mX m 

= i V • (5,Xu) + V • (XV^?) + ^ (a,,8) • B 

X 4mX m 

~V.[X(V^)sin^C] + ^(SxB). 

=^[Sxa,(Xa^8)], + ^(8xB)„ (A13) 

where we substituted Eq. (IA9E Considered together as a vector equation, Eqs. (IA10p - (IA13l) coincide with Eq. (l47)) . 
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